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In first-order inflation a phase transition is completed by the collisions 
of expanding true-vacuum bubbles. If bubble collisions produce large num- 
bers of soft scalar particles carrying quantum numbers associated with a 
spontaneously broken symmetry, then symmetry restoration may occur in 
a "pre-heating" phase in a manner similar to symmetry restoration in the 
pre-heating phase of slow-roll inflation. Since bubble collisions lead to inho- 
mogeneities, there is the possibility of inhomogeneous symmetry restoration 
where restoration occurs only in the regions of wall collisions. 
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To a very good approximation the universe was in local thermodynamic equilibrium 
(lte) for nearly all of its early development. However, there should have been brief, but 
important, departures from LTE. These excursions from equilibrium left an imprint on 
the universe. Examples of non-LTE phenomena include baryogenesis, nucleosynthesis, 
freeze-out of a massive particle species, decoupling of matter and radiation, production 
of topological or non-topological defects in cosmological phase transitions, inflation, and 
reheating after inflation. In fact, it may be argued that nearly all of early- universe cos- 
mology is the study of departures from lte. It is commonly believed that many of the 
current issues in cosmology require an understanding of the nontrivial dynamics in the 
approach to equilibrium in the early universe. Nevertheless, despite its immense rele- 
vance, only very recently has substantial effort been devoted to a detailed understanding 
of nonequilibrium phenomena in the early universe. 

The non-equilibrium process of interest in this study is the phenomenon of reheating 
after inflation. There are many varieties of inflation models, but all have an early period 
of rapid expansion of the universe where the Robertson- Walker scale factor 'accelerates' 
(i.e., a > 0). At the end of the accelerated-expansion phase the radiation density of the 
universe is effectively zero, and the universe must be 'reheated'. [] 

In 'slow-roll' (sometimes referred to as 'chaotic') inflation models Q, the universe 

after inflation was dominated by the energy density contained in the coherent motion 

of a scalar field known as the inflaton, whose potential energy density was responsible 

for the accelerated expansion. Reheating in slow-roll inflation involves conversion of this 

coherent scalar-field energy density into into a thermal distribution of radiation. In a 

simple scenario of reheating, the inflaton field coherently oscillated about the minimum 

of its potential until the age of the universe was equal to the lifetime of the inflaton, then 

1 Of course 're-heated' may be somewhat of a misnomer since there is no guarantee that the universe 
was hot before inflation. 
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the inflaton decayed, and the decay products thermalized. 

Recent investigations into the non-linear quantum dynamics of scalar fields have 
implications for reheating after slow- roll inflation [[!]]. These studies reveal that the 
scenario by which the energy density in coherent oscillations of the inflaton field is 
converted to radiation may differ significantly from the above picture, which considered 
only the linear evolution in time of the inflation field . Quantum nonlinear effects may 
lead to an extremely effective dissipational dynamics and explosive particle production 
in even the simplest self-interacting theory where single particle decay is kinematically 
forbidden. It is possible that almost all of the energy stored in the form of coherent 
inflaton oscillations at the end of inflation is released after only a few oscillation periods. 
The energy is released in the form of inflaton decay products, whose occupation number 
is extremely large, and have energies much smaller than the temperature that would 
have been obtained by an instantaneous conversion of the inflaton energy density into 
radiation. 

Since it requires several scattering times for the low-energy decay products to form 
a thermal distribution, it is rather reasonable to consider the period in which most of 
the energy density of the universe was in the form of the non-thermal quanta produced 
by inflaton decay as a separate cosmological era. This is generally referred to as the 
'preheating' epoch. 

The phenomenon of symmetry restoration during the preheating era has been in- 
vestigated recently by Tkachev || and by Kofman, Linde, and Starobinski Q in the 
framework of typical chaotic inflationary models. It was shown that symmetry restora- 
tion processes during the nonequilibrium stage of preheating may be very efficient with 
important implications for Grand Unified Theories (guts) and axions. Indeed, if a GUT 
symmetry is restored during the preheating epoch, the subsequent symmetry breaking 
phase transition will reintroduce the problems of monopoles || or domain walls [|]. 
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In first-order inflation models (generally, any model in which inflation is completed by 
a strongly first-order phase transition, e.g., the extended inflationary scenario proposed 
by La and Steinhardt [|7|]; for a review of first-order inflation models, see ||) the universe 
was dominated by scalar-field vacuum energy as in slow-roll inflation, but inflation was 
terminated by the nucleation of true vacuum bubbles. At the end of first-order inflation 
most of the energy density of the universe was contained in the bubble walls. Reheating 
was instigated by the collisions of bubble walls, which converted the bubble-wall tension 
into individual quanta of the scalar field, which then decayed into normal particles, which 
eventually scattered and formed a thermal distribution. 

The aim of the present paper is to suggest another situation in which symmetry 
restoration can occur efficiently out-of-equilibrium, namely during the preheating era 
subsequent to first-order inflation. 

As discussed above, the basic idea of reheating in first-order inflation is essentially 
the same as in chaotic inflation: energy initially stored in a coherent scalar field must 
be converted into radiation. However, in first-order inflation this releasing of energy 
takes place through a number of steps involving both classical and quantum processes, 
and a rich phenomenology associated with these scenarios can arise. For example, it 
has been suggested that gravitational waves [p|,[l0|1 , black holes HiTjp!2| and the baryon 



asymmetry |l3j may have been produced during the phase transition. Whether or not 



such phenomena actually occur depends in part on the details of reheating. For instance, 

in the baryogenesis scenario of Ref . |1| it is important to know if the only source of heavy 

GUT bosons is from primary particles produced in the bubble wall collisions which, in 

turn, depends crucially whether the GUT symmetry is restored after bubble collisions, 

i.e., on the value of the re-heat temperature, Trh-Q 

2 The reheating temperature, Trh, is usually defined as the temperature of the universe when the 
thermal spectrum of radiation was first obtained after inflation. 



We shall show, however, that similar to what occurs in the chaotic inflationary sce- 
narios, the details of symmetry restoration may turn out to be rather independent of 
Trh, and may in fact be quite complicated, with the symmetry restored in some regions 
of the universe, but not others. 

In order to keep the discussion as general as possible, we will not specify any particular 
first-order inflaton model, but describe the salient features of the inflaton potential in 
terms of three parameters (A CT , a , and e). We denote the inflaton field by a, which 
has a potential of the general form suitable to provide for a first-order phase transition. 
(Table 1 lists the fields and their interactions.) The potential will be described in terms 
of a dimensionless coupling constant X a , a dimensionless constant e that determines the 
splitting between false-vacuum and true-vacuum potential energy densities, and a mass 
scale o"o, which also plays the role of the vacuum expectation value when the symmetry is 
broken. The mass of the field will be \]/ 2 cro, and the difference in energy density between 
the false and true vacuum states will be denoted as AV = e\ a o-Q. The parameter e must 
be less than unity for sufficient inflation to occur. This also implies that the bubbles 
of true vacuum formed in the transition will be "thin-wall" bubbles, with wall thickness 
much smaller than the radius. 

From the few parameters X a , e, and a , one can find all the information required about 
the bubbles formed in the phase transition. For instance, in the thin-wall approximation, 
the size of a nucleated bubble is given by R c ~ {e\ 1 J 2 a^j . Bubbles with a radius 
smaller than this critical size will not grow, whereas bubbles larger than the critical 
size are exponentially disfavored. Another crucial parameter is the thickness of the 
wall separating the true-vacuum region inside from the false-vacuum region outside the 
bubble: A ~ (\)J 2 (jQj . The ratio of the bubble-wall thickness to its size is A/Rq ~ e, 
which is much less than unity if the thin-wall approximation is adopted. Finally, the 
energy per unit area of the bubble wall is r] ~ \ l J 2 al. 
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Table 1: Three fields are involved in our consideration: the inflaton field a; the field x into 
which the domain walls disperse; and 0, a field whose spontaneously broken symmetry 
may be restored by the x background. In some models x an d cr may be the same field. 



Interaction 


Potential term 


inflaton self interaction: 


V (a) 


= K(v 2 - 




inflaton false-vacuum energy density: 


AV = 






X — 4> interaction: 


V x <f> = 


g<t> 2 x 2 




(p self interaction: 


V (<t>) 


= h^ 2 - 





When a bubble wall forms, false-vacuum energy is transformed into bubble-wall en- 
ergy, with the wall energy initially in the form of static surface energy. As the bubbles 
expand converting false vacuum to true vacuum, more and more of the wall energy 
becomes kinetic as the walls become highly relativistic. Numerical simulations [|TQ|JTT|] 
demonstrate that during collisions the walls oscillate through each other, dispersing the 
kinetic energy at a rate determined by the frequency of these oscillations. When the 
bubbles have slowed after a few oscillations, they then dissipate their surface energy into 
particles of typical energy determined by the wall thickness. 

Although the particles produced in the initial collisions of the walls may play an 
interesting role in preheating and reheating, in the following we will concentrate on the 
implications of the particles produced by the potential energy density of the bubble walls. 
Bubble walls can be envisaged as coherent states of inflaton particles, so that the typical 
energy of the products of their decays is simply the mass of the inflaton. This energy 
scale is just equal to the inverse thickness of the wall. 
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Let's envision the collision of two plane-parallel domain walls. The potential energy 
per unit area of the bubble walls is given by 77 ~ Ay 2 o"o. Taking the mean energy of the 
particles produced in the bubble wall collisions to be of order of the inverse thickness of 
the wall, E ~ A -1 , the mean number-per-area of particles produced from the potential 
energy in the collisions is N ~ rj/E ~ \)J 2 <JqA. 

Let's now assume that the particles are spread out a distance d from the region of 
the wall collision. If we approximate the particle density as uniform out to a distance d, 
then the particle number density within the region is simply 

n = N/d~ Xl /2 alA/d ~ aljd. (1) 

In the limit that the walls are spherical with radius R and the collision products instantly 
fill the bubble interior, then the factor of d in Eq. (|I|) should be replaced by R. 

Eventually the products of bubble-wall collisions will be redistributed throughout 
the bubble interior and thermalized. If we assume that thermalization is instantaneous, 
the reheating temperature is found by imposing pn = ((7*7r 2 /30)T^ H = Ay, where is 
the effective number of degrees of freedom in all the species of particles formed in the 
thermalization processes. Using AV = eAo-cx 4 results in a re-heat temperature of Trh ~ 
gl 1 ' 4 e 1 / 4 \ l J A a . Let us now assume that the typical energy of the particles produced 
through bubble collisions is smaller than Trh, i-e., A -1 < Trh, which translates into the 
condition (taking g* ~ 100) X a < 10~ x e. If this condition is satisfied, then a period is 
required for equilibration, namely for particles to scatter from energies approximately 
equal to A -1 to a thermal distribution of temperature Trh- In addition, since the bubbles 
were originally empty, homogenation is not instantaneous, and requires a time at least 
as long as the light travel time across a bubble. If either of these two time scales is 
sufficiently long, we may consider the time interval during which particles do not have a 
homogeneous thermal distribution function as a separate epoch: the preheating era. 
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As a first approximation, during the preheating period the distribution function of 
the created particles can be chosen of the form |3j 



where E = A -1 and the constant A may be fixed by computing the number density of 
particles, n = (2ir)~ 3 f cPp f(p), and setting it equal to the estimate given in Eq. ([[]). Of 
course, A has mass dimension one. 

Let us now imagine that particles x are produced in the bubble wall collisions and are 
charged under some symmetry group, so that their mass m x depends upon some scalar 
field <p as m 2 (<fi) = rriQ + g<p 2 -^ Here, g represents a combination of numerical factors and 
a coupling constant. As a simple example we might assume that the ^-dependent mass 
originates from a potential term of the form V X( f, = g<f) 2 x 2 - 

As opposed to large-angle scattering processes, forward-scattering processes do not 
alter the distribution function of the particles traversing a gas of quanta, but simply 
modify the dispersion relation. This remains true also in the case of a nonequilibrium 
system. Forward scattering is manifest, for example, as ensemble and scalar background 
corrections to the particle masses. Since the forward scattering rate is usually larger than 
the large-angle scattering rate responsible for establishing a thermal distribution, the 
nonequilibrium ensemble and scalar background corrections are present even before the 
initial distribution function, Eq. (0), relaxes to its thermal value. These considerations 
allow us to impose uj 2 = p 2 + m' x ((j)) as the dispersion relation for the particles created 
by bubble collisions. 

We can not use the imaginary-time formalism to determine the effective potential for 
the scalar field during the nonequilibrium preheating period since in the nonequilib- 
rium case there is no relation between the density matrix of the system and the time 
3 Of course by <p 2 and \ 2 we mean the appropriate sum over the members of the group representation. 




(2) 
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evolution operator, which is of essential importance in the formalism. There is, how- 



ever, the real-time formalism of Thermo Field Dynamics, which suites our purposes [|T4 |. 
The contribution of the particles created by bubble collisions to the one-loop effective 
potential of the scalar field can be written as 

The first integration in u must be done treating u as a free parameter and setting 



u! p (4>) = J p 2 + m 2 (0). By making use of Eq. (g), one obtains 

AV(0) = -A J -0^ 6[E - u, p (0)] ~ I [ml (0) - £ 2 ] . (4) 

Since we are interested in the 0-dependent part of the potential, we can ignore the 
nE term and the factor of m\ in m 2 ((f)), and write the potential for the non-equilibrium 
configuration as AV(<p) = B NE (fi 2 , where B NE = gn/E. A similar expression was obtained 
by Tkachev in Ref. ||, using the definition of the effective potential as (the negative of) 
the pressure of the system, and assuming that the number of particles does not change 
on time scales of interest as the field evolves. 

We now use n = o^jd from Eq. ([I]), and E ~ A -1 , to obtain B NE ~ ga^A/d. Of 
course d will depend upon the details of the model and the complexities involved in the 
completion of the phase transition. But it is reasonable to expect, at least initially, that 
d is of order A, so let us write d = £A. Of course as the bosons diffuse into the bubble 
interior £ will change in time, so we expect £ to grow and eventually to become much 
greater than unity. But initially, at least, £ should not be too much larger than unity. 
In terms of £, we may express _E?ne as B NE ~ g<?o/£,- 

Now there are two things left to do. First, we will determine the conditions un- 
der which the non-equilibrium contributions to the effective potential can restore the 
symmetry, and then determine the criterion for the non-equilibrium effects to be more 
important than the equilibrium effects obtained after re-heating. 
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Let us take the <fi tree-level self-interaction potential to be of the form Vo(0) = 
<V(</> 2 — <Po) 2 - The symmetry will be restored (i.e., = will be a stable minimum) if 
d 2 V/dcf) 2 evaluated at </> = is positive, where now V includes the sum of the tree-level 
potential and the one-loop correction, V = Vq + AV. Symmetry restoration will occur 
due to non-equilibrium effects if — A^ 2 , + -^ne > 0. This translates into a bound on £ 
for symmetry restoration: 

We can imagine three interesting limits depending upon the magnitude of the left-hand 
side (lhs) of this inequality. Since we expect £ always to be greater than one, if the 
LHS is less than unity we would expect non-equilibrium effects never to cause symmetry 
restoration. If the lhs is greater than one but not very large, then one might expect 
temporary restoration of symmetry around the regions of bubble collisions. Then as £ 
starts to grow as the bubble interior is filled, the symmetry will be broken when the 
inequality is violated. Finally, the lhs may be so much greater than unity that the 
symmetry is restored even after the bubble interiors are filled. 

Of course the symmetry may also be broken after re-heating by thermal effects. 
This can be seen by calculating the ^-dependent term in the one-loop effective potential 
obtained by assuming that the system is in lte at temperature Trh- Including V x ^ = 
g<p 2 X 2 , in the high-temperature limit the one-loop thermal corrections lead to AV((f), T) ~ 
gT 2 <)) 2 + A TV- If we write AV(<f>, T RH ) = B EQ <p 2 , with B EQ = (g + A^)T RH 2 , then 5 NE 
plays a role in non-equilibrium transitions similar to that played by I?eq for thermal 
transit ions. Q 

4 Thus, we see that so far as symmetry restoration is concerned, in the presence of the soft bosons left 
behind in the debris of wall collisions, a scalar held behaves as if it was in lte at an effective temperature 
T£ FF - B NE /(g + A ) ~ alg/[i{g + A )]. 
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Symmetry will be restored after re-heating if —X^q + -Beq > 0, or expressing this as 
a limit to Trh: Trh 2 > X<f>4>o/(9 + -V). Now we know Trh in terms of the parameters 
of the inflaton potential, so we may express the criterion for symmetry restoration after 
re-heating as 

A, $ V eK ' { ) 

The condition for symmetry restoration in pre-heating, Eq. (El), and the condition 



for symmetry restoration in re- heating, Eq. (f|), are most easily contrasted in the limit 
g > \<f,. In that limit 

9 f £ (symmetry restoration during pre-heating) 

f S> , (■) 

* I ^Jg*/e\ a (symmetry restoration during re-heating). 

Depending upon the parameters, it is possible to have restoration during both pre- 
heating and re-heating, during neither pre-heating or re-heating, or during one and 
not the other. Of particular interest might be the case where restoration occurs only 
during pre-heating when £ is not too large. Then the effects of inhomogeneous symmetry 
restoration will not be erased during re-heating. 

In conclusion, symmetry restoration may well occur in the preheating phase following 
first-order inflation. Unlike symmetry restoration in the preheating phase of chaotic 
inflation, the restoration may be inhomogeneous after first-order inflation. The basic 
point is that the phase-space density of bosons created in wall collisions is greatest in 
regions of wall interactions. One may imagine situations where restoration occurs among 
the debris of wall collisions, but not in the initially empty interior of the bubbles. In 
such a case, the subsequent symmetry breaking restoration might result in creation of 
topological defects if the region of wall interactions is large enough to contain these 
defects. 

Cosmological implications of this possibility require further study. 
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